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We present a detailed theoretical study of the ultrafast quasiparticle relaxation dynamics observed 
in normal metals and heavy fermion materials with femtosecond time-resolved optical pump-probe 
spectroscopy. For normal metals, a nonthermal electron distribution gives rise to a temperature (T) 
independent electron-phonon relaxation time at low temperatures, in contrast to the r~^-divergent 
behavior predicted by the two-temperature model. For heavy fermion compounds, we find that the 
blocking of electron-phonon scattering for heavy electrons within the density-of-states peak near the 
Fermi energy is crucial to explain the rapid increase of the electron-phonon relaxation time below 
the Kondo temperature. We propose the hypothesis that the slower Fermi velocity compared to the 
sound velocity provides a natural blocking mechanism due to energy and momentum conservation 
laws. 
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I. INTRODUCTION 

Recently, ultrafast time-resolved optical spectroscopy 
has been used to reveal the nature of the quasipar- 
ticle relaxation dynamics in condensed matter sys- 
tems. Femtosecond time-resolved pump-probe opti- 
cal measurements have been carried out in normal 
metalsfi*^ conventional^ and high-Tc superconductors;* 
and charge-density wave solidsi^ Time-resolved terahertz 
spectroscopy has measured picosecond time-scale tran- 
sient conductivity in colossal magnetoresistance mangan- 
itesi^ These experiments show directly in the time do- 
main how constituent degrees of freedom of materials 
interact with each other, which is important to under- 
stand the physics governing the ground state and the 
low-energy excited states of materials. Understanding 
the fast dynamics of quasiparticles is also crucial for tech- 
nological applications of these materials, for example, in 
the design of very fast switching devices. 

We recently reported femtosecond time-resolved 
pump-probe optical measurements on LuAgCu4 and 
YbAgCu4.- These two materials are isostructural with 
a negligible difference in lattice constants* (about 0.3 %) 
and about 1 % difference in atomic mass between Lu and 
Yb. However, their electronic structures are strikingly 
different due to the /-level occupancies in Lu and Yb. 
The closed-shell 4/ levels in Lu have no significant inter- 
action with Cu 3d and Ag Ad conduction electrons, and 
LuAgCu4 is a normal metalA In contrast, the open-shell 
4/ levels in Yb, that is, the localized hole in the 4/^^ 
configuration, have a strong interaction with conduction 
electrons, and make YbAgCu4 a heavy- fermion material, 
which is characterized by a large peak in the electron den- 
sity of states (DOS) at the Fermi energy (Ep) and, equiv- 
alently, a large Sommerfeld coefficient 7ii2iii Our time 
resolved optical experiments show very different relax- 
ation dynamics in these two materials. LuAgCu4 shows 
a relaxation time (r) versus temperature (T) behavior 
similar to other normal metals, such as Ag and Au.i^ In 



particular, there is very little T dependence in r at low 
temperatures. However, YbAgCu4 shows approximately 
a 100-fold increase in r, as T is decreased from the Kondo 
temperature (Tk) down to 10 K. Kondo temperature is 
about 100 K in YbAgCu4 and typically characterizes the 
width of the large DOS near Ep. 

In Ref. 0, we presented experimental data showing dif- 
ferent relaxation dynamics in LuAgCu4 and YbAgCu4 
along with the main theoretical ideas and final results of 
the calculations based on coupled electron and phonon 
Boltzmann transport equations. In this paper, we re- 
port on the details of the theoretical model and analysis, 
and clarify the underlying physics. The theories that 
existed before our study are summarized in Sec. II A. 
Sections II B and HI discuss our calculations of the re- 
laxation dynamics in normal metals and heavy fermion 
materials, respectively. Section IV summarizes our re- 
sults. 



II. NORMAL METALS 

A. Two-temperature model, its inconsistency with 
experiments, and the nonthermal electron model 

In ultrafast optical pump-probe spectroscopy, an ul- 
trafast laser pulse initially excites the electron system, 
and the probe pulse monitors the relaxation of the elec- 
tron system by measuring transient optical properties 
with subpicosecond time-resolution. Because the diffu- 
sion of heat out of the illuminated (probed) region is 
much slower (tens to hundreds of nanoseconds) than the 
time scales of interest, the relaxation of the excited elec- 
trons is due to the thermalization among electrons and 
other degrees of freedom, such as phonons, within the 
illuminated areas. 

The most commonly used model for the relaxation dy- 
namics of photoexcited electrons in metals is called the 
two-temperature model (TTM)jiS which assumes much 
faster relaxation within each constituent system (e.g., 
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electron system, phonon system) compared to the relax- 
ation between these constituent systems. In this approx- 
imation, the temperature of each system can be defined 
during relaxation, and the relaxation time r between sys- 
tem 1 and 2 is determined by their specific heats, Ci and 
C2, and the energy transfer rate per temperature differ- 
ence, g(T)Mi 
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Kaganov, Lifshitz, and Tanataro\ii^ calculated the en- 
ergy transfer rate g{T) between electrons (e) and phonons 
(p) in normal metals by solving coupled Boltzmann trans- 
port equations for electrons and phonons with ther- 
mal equilibrium distributions at different temperatures. 
Their results predict t ^ T dXT > Tjjebye/^ and r ~ 
at T < Tuebye/^, whcrc Tuebye IS the Debyc temperature. 

Groeneveld et almi^ have measured the relaxation time 
T in Au and Ag as a function of temperature (T = 
10 — 300 K) and laser-energy density {Uiaser = 0.3 — 1.3 
J cm~'^) using the femtosecond pump-probe technique. 
Their experimental data was inconsistent with the TTM 
predictions for r versus T and r versus Uiaser, which led 
Groeneveld et al. to consider the nonthermal electron 
model (NEM). Using a simple estimate, they pointed 
out that the relaxation time within the electron system 
is comparable to the electron-phonon relaxation time. 
Instead of assuming a thermal equilibrium distribution 
for the electrons, they numerically solved the Boltz- 
mann equation for electrons with electron-electron and 
electron-phonon scattering, starting from an initial non- 
thermal electron distribution created by a laser pulse. 
(Based on the fact that the phonon specific heat is much 
larger than the electron specific heat, phonons were as- 
sumed to have a thermal distribution with a time inde- 
pendent T.) This simulation could explain experimen- 
tal data down to about 50 K. Although Groeneveld et 
aL's work proposed the essential idea of a nonthermal 
electron system, their analysis focused on temperatures 
above 50 K, and excluded the low-T region where the 
most striking difference between the TTM (t ~ T~^) 
and the experimental data (almost T-independent r) oc- 
curs. Therefore, we undertake a more detailed analysis, 
particularly focused on the low-T region, to understand 
the difference in relaxation dynamics between thermal 
and nonthermal electrons in normal metals as well as in 
heavy fermion metals. 



B. Relaxation dynamics between nonthermal 
electrons and phonons 

We consider the coupled Boltzmann equations for 
both electrons and phonons, with electron-electron and 
electron-phonon scattering included»i^ The Boltzmann 
equations with momentum indices are transformed into 
the following equations with energy indices for a model 
with isotropic Debye phonons and electrons with an 



isotropic parabolic dispersion relation. In terms of elec- 
tron and phonon distributions at time t, fe(t) and b,^(t), 
and the electron and phonon DOS, De{e) and Dp{ijj), 
where e and uj represent the electron and phonon ener- 
gies, the coupled Boltzmann equations araiSiiSiii 
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The linear dependence of the electron-phonon scatter- 
ing rate on the acoustic phonon energy uj [see, e.g., 
Eq. (26.42) in Ref. has been canceled out by a l/w 
factor which originates from changing momentum indices 
to energy indices, so that K^p is an energy independent 
constant. Kep and K^^ include the squares of the scat- 
tering matrix elements and all other numerical factors, 
and have the units of [energy/time]. 

In the remainder of this section, we consider only light 
electron systems (normal metals) . The above Boltzmann 
equations, with electron and phonon DOS explicitly in- 
cluded, show why the electron-electron relaxation time is 
comparable to the electron-phonon relaxation time at low 
temperatures. The relaxation rate depends on the square 
of the scattering matrix elements (i.e., Kf,p and -fCee) and 
the number of available final states (i.e., De and Dp). If 
we compare Eqs. @ and ©, even though Kf.e is typically 
much larger than Kep, the large disparity between the 
magnitude of the phonon {Dp) and electron [D^) DOS 
expedites the relaxation of electrons with phonons, and 
hinders electrons from reaching thermal equilibrium. The 
electron-electron and electron-phonon relaxation times. 
Tee and Tep, vary with T. In Ref. ,2, Groeneveld et al. 
found Tee ~ T'^ [Eq. (16) of Ref. from Fermi hquid 
theory. Figurenschematically compares the temperature 
dependence of Tee (thick line) and Tep (thin line) , the lat- 
ter of which is from the TTM results. ^'^^ The plot shows 
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FIG. 1: Schematic plot for the temperature dependence of 
electron-electron and electron-phonon relaxation times, Tee 
(thick line) and rep (thin line). Tee is from Eq. (16) in 
Ref. Tep is from the TTM calculation, which assumes a 
thermal electron distribution. The ratio between Tee and rep 
at high T (e.g., Toetye or 300 K) depends on the Kee/Kep 
ratio. For Kse/ Kep=700 used in our calculation, Tee < Tep at 

T = Toebye: ^ shoWU iu the plot. 

that it is likely that at high temperatures the electron 
system thermalizes within itself faster than the thermal- 
ization time with the lattice, whereas at low temperatm'es 
electron-electron thermalization becomes slower than the 
electron-phonon relaxation. More detailed discussion on 
T-dependent thermal versus nonthermal electron distri- 
butions will be given later in this section. 

We now discuss the methods and parameter values we 
used to analyze the above Boltzmann equations. The 
coupled Boltzmann equations are solved in two different 
ways. In the first method, an initial condition at t — is 
set as a nonthermal electron distribution excited by the 
laser pulse, ff{t = 0) = fpoif^^Ti) + Af^, and a thermal 
phonon distribution, buj{t = 0) = bBEiy>,Ti), where T^, 
fFD{^,Ti), and bBE{^^,Ti) represent the initial temper- 
ature, Fermi-Dirac, and Bose-Einstein distribution func- 
tions before the photoexcitation. We find f^{t) and b^{t) 
at subsequent times by solving the differential equations 
until both /c (t) and bai (t) approach the final thermal dis- 
tributions at equal T's. A small A/^ has been chosen 
so that the final T (Tf) is less than 3 K higher than Ti 
(even at the lowest T), which simulates the low laser in- 
tensity used in the experiments^ In the simulation the 
total electron energy Ee{t) is calculated at each time step, 
and the instantaneous relaxation time rit) is evaluated 
using Ee{t) at three consecutive time steps. 

In the second method, we perform a linear stabil- 
ity analysis around the final states, /F_D(e,T^) and 
6s_e(w,T/), by expanding the coupled Boltzmann equa- 
tions, Eqs. ||2J) and (|3J), linearly in Sfe{t) and 6buj{t), 
where /,(i) = fpoie^Tf) + Sf,{t) and b^t) = 
bBE{i^,Tf) + 6buj{t). The zeroth order terms of the ex- 
pansion vanish on both sides of the equations, because 
they correspond to the final equilibrium state. With 
Sf^{t) = ufe"*/^, Sb^{t) = yPe-*/^, and discretized N 
electron and M phonon energy levels, the linear differ- 



ential equations can be cast into an eigenvalue problem 
of an {N + M) x {N + M) matrix, the solution of which 
gives — l/r (eigenvalues) and the normal modes of the 
relaxation. Two of the normal modes have an unphysi- 
cal infinite relaxation time r, which originates from to- 
tal energy and total electron number conservation. The 
rest of the modes represent all possible relaxation modes 
of the system. In general, the relaxation of the sys- 
tem can be represented as a linear combination of these 
modes. By examine the eigenvector of each mode, we 
find whether the mode predominantly contributes to e- 

e, p-p, or e-p relaxation. Since the energy transfer rate 
from the electron system to the phonon system at time 
t of a specific mode with relaxation time r is given by 
dEe,rit)/dt = [/(e - EF)i-l/T)vtDe{e)de]e-*^'' , a use- 
ful quantity to identify the e-p relaxation modes is the 
following e-p energy transfer strength rsii)'- 

1 ^ 

rE{i) = ^ - i5i.)<(i)De(e„)Ae, (7) 

where Ae is the energy step size, and 
(wf , ^2, f^, f 1 , is the normalized eigenvector 
for mode i with proper overall sign. r^(i) characterizes 
the effectiveness of mode i for electron-phonon energy 
relaxation. The modes which have large te dominate 
in e-p relaxation, and their eigenvectors describe how 
the relaxation between the electron and phonon systems 
occurs. The modes with small ^^(i) are either p-p or 
e-e relaxation modes, and describe how the relaxation 
within each system happens. 

The electron and phonon energy levels are discretized 
with a step size of 2 meV, which we find small enough 
that any smaller step size would not change our results 
even at 10 K (the lowest T of the calculations and exper- 
iments). The energy window for the electrons has been 
chosen to be between -0.15 eV and 0.15 eV with Ep = 0, 
which is wide enough in comparison to the highest T 
of the calculations ('^ 300 K). Since the typical band 
width of normal metals is of the order of eV, the elec- 
tron DOS for normal metals is assumed to be constant 
within this energy window. Fitting the experimental elec- 
tronic specific heat (Ce) data, we obtain — 2.1 eV^^ 

f. u.^^ spin~^ for LuAgCu4jSi which has been used for all 
of the results presented in this section. Phonons are mod- 
eled by the Debye phonon model. Since only longitudinal 
phonon modes couple with electrons in the isotropic elec- 
tron and phonon modelji^ we use the longitudinal phonon 
DOS with a Debye energy ujd — 24 meV (Ref. Il0|) and 
6 atoms per unit cell, Dpiuj) — 18a;^/w|), in the Boltz- 
mann equations. For most of the calculations, we use 
Kee/Kf,p = 700 (see Figs. 121121 and related discussion for 
this choice of Kee/Kep ratio), and K^p has been chosen 
as Kep — 0.93 eV/ps for normal metals by fitting ex- 
perimental data with the results of the calculation (see 
Fig.EJ. 

We now discuss the results obtained from the first 
method. Figure |21 shows a typical result of normalized 
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FIG. 3: Relaxation time r calculated from the coupled Boltz- 
mann equations for LuAgCu4, along with experimental data 
and the TTM prediction.^ 



FIG. 2: A typical result of (a) normalized excess electron 
energy versus time and (b) the instantaneous relaxation time 
versus time at high T (solid lines, T = 18 meV = 208 K) and 
low T (dotted lines, T = 4 meV = 46 K). 

excess electron energy, ASg/AiJe.maa; versus time (i), 
and r (instantaneous relaxation time) versus t at high 
(sohd Hues) and low (dotted hues) temperatures, where 
A£;e(t) = E^it) -Ee{t = oo) and /^E^.^ax = AEe{t = 
0). Two different stages of relaxation can be identified 
for high T results. During the intermediate time be- 
tween 0.3 ps and 2 ps (as mentioned earlier, the timescale 
in picosecond is obtained by using K^e/Kep = 700 and 
i^ep — 0.93 eV/ps), the energy transfer from electrons 
to phonons is fast. After more than 99 percent of the 
excess electron energy is transfered to phonons, a slow 
relaxation (about 100 times slower than the intermediate 
stage) appears. The changes of distribution functions at 
different time stages show that the fast relaxation dur- 
ing the intermediate time corresponds to electron-phonon 
relaxation, whereas the later slow relaxation corresponds 
to phonon-phonon thermalization processesJ^ The relax- 
ation behavior is independent of the initial conditions, 
i.e., the form of A/^, except at very early stages of the 
relaxation, namely, t < 0.3 ps in Fig.|2| Since the changes 
in electron distribution during the late stage are expected 
to be unobservable due to experimental noise, we identify 
the intermediate stage t [t = 0.3 ~ 2 ps for the solid line 
in Fig.[5fb)] as the one measured experimentally at high 
T. As T is lowered in the simulation, the well-defined flat 
intermediate time region in the t versus t plot is replaced 
by a gradually changing r with a minimum as shown in 
Fig.|3b). This indicates multiple relaxation times, a sign 
of a nonthermal electron distribution. As done in the ex- 
periments, we find the best fitting single relaxation time 
through a linear fit of log \dEe/dt\ versus t, which is close 
to the minimum t of the dotted line in Fig.[2Ib). 

The r's obtained at various T's are plotted in Fig. 
(solid circles) along with the experimental data (open 
circles) for LuAgCu4 and the TTM prediction (dashed 



line, see Ref. The time unit has been scaled to fit 
the experimental data at high T, which gives K^p — 
0.93 eV/ps as mentioned previously. The fit does not 
uniquely determine K^e, as long as K^e is smaller than 
or similar to 700 x Kf.p- The results show reasonable 
agreement with experimental data, including the region 
below 50 K, where the TTM predicts strikingly different 
T - T-^ behavior (Refs.H and[l2). By comparing df / dt 
and —d/po/dT, we examine whether the electron system 
approaches the final equilibrium state while maintaining 
a thermal distribution or not, as shown in Fig. ^ at high 
and low T's. At high T [Fig. ^a)], the two curves co- 
incide with each other [indistinguishable in Fig. Ela)], 
implying that the electron system has a thermal distri- 
bution. However, at low T [Fig. ^b)], df/dt (solid line 
with solid circles) has a width of the order of the Debye 
temperature (Tuebye) instead of T, which is the width of 
—dfro/dT (dotted line)i2£ The results clearly show that 
at low temperatures the electron system does not have 
a thermal distribution, which agrees with the previous 
discussion in relation to Fig. ^ 

The effect of a nonthermal electron distribution on the 
relaxation dynamics becomes even more evident, if we 
increase K^e to 7000 x Kep so that the electron-electron 
relaxation is faster than the electron-phonon relaxation 
even in the low-T region. The result (solid circles) is 
shown in Fig. along with TTM prediction (line) . In 
this case, the electron system has a thermal distribution 
in the whole T range, and the TTM prediction r ^ T^'^ 
is recovered at low T. It is remarkable that our simu- 
lation with a large K^e recovers the TTM prediction in 
spite of the fact that our Boltzmann equations do not 
include direct phonon-phonon scattering and, therefore, 
the phonon distribution is nonthermal except at i = oo. 
We speculate that the very weak dependence of e-p ther- 
malization dynamics on the phonon distribution (in con- 
trast to the very strong dependence on electron distribu- 
tion) originates from the characteristics of Bose (versus 
Fermi) statisticsi^ 
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FIG. 4: Normalized df/dt compared with normalized 
-dfpD/dT at (a) high T and (b) low T. For high T, the 
two lines coincide. 



We provide a qualitative reason why the thermal elec- 
tron distribution gives rise to a slow r T^^ relaxation 
behavior and the nonthermal electron distribution gives 
faster and less-T-dependent relaxation behavior at low T . 
As depicted in Fig. EJa) , if the electrons have a thermal 
distribution at T (solid line), which is slightly higher than 
the final T (dotted line), the electron-phonon scattering 
important for the relaxation happens within the energy 
range of the order of T from Ep. Therefore, the relax- 
ation rate depends on how many phonon modes exist at 
u! < T. In the Debye phonon model, since Dp{uj) ~ w^, 
the relaxation rate is proportional to T^. In con- 
trast, if the electron distribution is nonthermal [solid line 
in Fig. |Sfb)], then electron-phonon relaxation occurs in 
over an energy range of the order of Debye energy [see 
also Fig. 0fb)]. This makes the electron-phonon relax- 
ation faster and less T-dependent, as indeed observed 
(Fig.©. 

In the second method, we linearize the coupled Boltz- 
mann equations and calculate eigenvalues and eigenvec- 
tors of the matrix. This second method supports and 
clarifies the results obtained earlier by the first method. 
With the same discretized energy levels as in the first 
method, a total of 163 energy levels (151 electron and 12 
phonon energy levels) exist in our model, which means 
a matrix of size 163 x 163 has to be numerically di- 
agonalized. Two of the 163 eigenmodes are unphysical 
as discussed earlier. The remaining 161 physical modes 
have negative eigenvalues, — 1/t, as expected for a stable 
fixed point, and represent all possible relaxation modes 
of the system. The r's at high and low T's are plot- 
ted in Figs. Ufa) andlJfc), respectively. The indices of 
modes i are assigned in descending order of r. The e-p 
energy transfer strength rp calculated by Eq. Q for high 
and low T's are plotted in Figs.^b) and[7|^d), showing 



Our model with K JK^p=7000 



TIM prediction 



10 H ^-3 ^ 

T ~ T at low T 



T(ps) 



10" 





10" 



100 200 
T(K) 



300 



FIG. 5; Solid circles: calculated relaxation time r from our 
model for LuAgCu4 for a very large K^e > 7000 x Kep- Line: 
TTM prediction from Fig. |3 With fast electron-electron re- 
laxation our simulations based on coupled Boltzmann equa- 
tions recover the TTM prediction r ~ T^'^ at low T. 
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FIG. 6: Schematic pictures explaining different relaxation 
dynamics for (a) thermal and (b) non-thermal electron distri- 
butions at low temperatures. Dotted lines represent the final 
thermal electron distribution. 
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the efficiency of eacfi mode for electron-pfionon energy 
relaxation. At both T's, there exist 11 modes (one less 
than the number of phonon energy levels) which have a 
much larger r than the rest of the modes. For example, 
at high T these modes are 2 orders of magnitude slower 
than the rest of the modes. These slow modes, whose 
eigenvectors have mostly phonon components and very 
little electron components, correspond to the late time 
phonon- phonon relaxation (see Fig. 13) . For the high T 
results in Figs.[7{a) and[7fb), a single mode with i = 12 
has a large and is very effective in transferring energy 
from the electron system to the phonon system in com- 
parison to the other modes. This mode has a value of r 
identical to the intermediate time t in Fig.|21found by the 
first method, and is well separated from other r values 
as seen in Fig. [Tfa).^^ The electron part of its eigenvec- 
tor matches —dfFD/dT, confirming the results obtained 
by the first method at high T [Fig. El^a)], that is, the 
presence of a thermal electron distribution. Modes with 
j=14, 41, and 158 have smaller te values than the i—12 
mode. Their eigenvectors show that they participate in e- 
p relaxation, but the electronic parts of the eigenvectors 
do not match —d/pD/dT. Their nonthermal behavior 
seems suppressed in the experiments and in the simula- 
tion due to the thermal i=12 mode, which has more than 
5 times larger energy transfer strength rg. The remain- 
der of modes between i=13 ^ 161 with very small te val- 
ues are mainly e-e thcrmalization modes. These features 
change as T is lowered, as shown in Figs.[7fc) and[7Jd). 
Many modes (modes with roughly te > 0.1 meV/ps) 
transfer energy from the electrons to the phonons. How- 
ever, none of their corresponding electronic eigenvectors 
can be described by —dfEo/dT, indicating the absence 
of a well defined electron temperature during relaxation. 
This again supports the results of the first method, pre- 
sented in Fig. 2fb). The rest of the modes between i—li 

161 with te very close to zero again correspond to e-e 
relaxation modes. Importantly, these have larger r val- 
ues compared to the high T case, showing the slowing 
down of e-e relaxation at low T, consistent with Fig. ^ 

Both methods show the importance of a nonthermal 
electron distribution for the electron-phonon relaxation 
dynamics at low T for normal metals. In the next sec- 
tion, we discuss how the relaxation dynamics of heavy 
fermion compounds are dramatically different from that 
of normal metals. 



III. HEAVY FERMION MATERIALS 

As mentioned in Sec. I, heavy fermion materials are 
characterized by a large DOS peak with a width of or- 
der Tk near Ep. YbAgCu4 has Tk ~ 100 K and is 
a paramagnetic metal down to the lowest T studied.-^'' 
The interaction between localized /-electrons and delo- 
calized conduction electrons in heavy fermion materials 
is a focus of intensive research in strongly correlated elec- 
tron systems and has not yet been fully understood. In 
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FIG. 7: [(a), (c)] Relaxation time r and [(b), (d)] e-p energy 
transfer strength re at [(a), (b)] high T and [(c), (d)] low T 
calculated by the second method, which linearizes the coupled 
Boltzmann equations. The index i is assigned in descending 
order of r. 



this paper, we consider the simplest model that captures 
the basic physics of heavy fermions, that is, a peak in 
the DOS near Ep- Even though some theories^^ pre- 
dict a disappearance of the peak DOS above Tk, we find 
that the calculated specific heat [Fig. |Htb)] and relax- 
ation time fFig. lll|l have only a weak dependence on the 
presence or absence of the peak DOS above Tk- (There- 
fore, neither the relaxation time data nor the specific heat 
data can distinguish whether the peak DOS is present or 
not above Tk ) In the hybridization gap model, local 
/ levels hybridize with the conduction band and open a 
gap with DOS peaks above and below the gap. If Ep 
is located within the peak, not in the gap, the results 
we will show below have little dependence on whether 
we use a hybridization gap model or a single peak DOS 
model. Therefore, in our simple model, we assume a 
T-independent single peak electron DOS as shown in 
Fig. ISJa). We further simplify the problem by choosing 
the Ep at the center of the peak, so that the chemical po- 
tential is T-independcnt. The peak DOS is described by 
a Gaussian function with a constant background DOS: 



De{e) ^ Dpeak exp[-(e/A)2] + Do 



(8) 



The discretized energy step size and energy window are 
identical to the normal metal case in Sec. II. We deter- 
mine the DOS parameters by fitting the electronic spe- 
cific heat data, Cg, as shown in Fig. |SIb), and obtain 
Dpeak = 70 eV^i f.u.-i spin-i, Dq = 2.1 eV-^ f.u.-i 
spin~^, and A = 13 meV. Dp^ak is directly related to the 
linear slope of Cg near T — 0. A is of the order of Tk 
and determines the temperature of the peak in Cg. Dq is 
identical to the electron DOS at Ep for LuAgCu4, and is 
related to Ce above Tk- Since LuAgCu4 and YbAgCu4 
are isostructural with almost identical atomic masses, we 
use identical phonon DOS for LuAgCu4 and YbAgCu4. 
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FIG. 8; (a) Model DOS and (b) calculated electron specific 
heat (solid circles) along with experimental data of Ce (line) 
for YbAgCu4. 



10-= 



10^ 



10' 

T(PS) 

10° 



10" 



Experiment 

Model with peak DOS only 



o o o oOqO : 



T ~ T at low T 



To gain insight into the nature of the scattering be- 
tween heavy electrons and phonons, we consider a case 
in which only heavy electrons with an isotropic parabolic 
dispersion exist. As in the normal metal case, we trans- 
form the coupled Boltzmann equations with momentum 
indices into the equations with energy indices. The re- 
sult is similar to the normal metal case, Eqs. Q-linjl, 
except for one important difference: If the Fermi ve- 
locity vp is slower than the sound velocity Vg, where 
vp = {de/dk)g^EF Vg = [duj / dk)k^Q, then the 

phonon integration in the electron-phonon scattering has 
a lower bound, which represents the blocking of electron- 
phonon scattering for low energy phonons and has a great 
influence on electron-phonon relaxation at low T . There- 
fore, if vp > Eqs. 0-10 with a peak DOS for Df, can 
be used to approximately model the relaxation dynamics 
for both light and heavy electrons, whereas if vp < Vs, 
appropriate blocking of scattering processes should be 
imposed upon these equations. 

We first consider the case vp > Vs, for which the im- 
portant physics is simply the increased DOS near Ep 
given by Eq. We use the same value of K^p and 

Kee as the LuAgCu4 case in Sec. II. The calculated 
T (solid squares) and the experimental data (open cir- 
cles) shown in Fig. [^disagree in two respects. First, the 
calculated r at around T=300 K is about 60 times less 
than the observed r. (Note that, experimentally, both 
YbAgCu4 and LuAgCu4 have similar r at ~ 300 K.) 
This difference is due to the increased electron DOS in 
the calculation, which enhances the e-p relaxation. An 
approximately 60 times smaller K^p would shift the cal- 
culated r in the whole T range by 60 times, but it is 
unlikely that K^p in YbAgCu4 would be smaller than 
Kep in LuAgCu4 in such a drastic way. Secondly, and 
more importantly, the divergence of calculated r's at low 
T is much weaker than the experimental data. This 
slow divergence in calculated r can be understood in 
the following way. The large electron DOS increases 
both the electron-electron and the electron-phonon scat- 
tering rates due to the increased number of available fi- 
nal states [see Eqs. Q-®]- However, electron-phonon 
scattering in Eq. Q has only one factor of D^, whereas 
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FIG. 9: Relaxation time calculated with the peak DOS only 
[Eq. Q and Fig. |HIa)], along with experimental data for 
YbAgCu4. 



electron-electron scattering in Eq. ijHJl has three factors of 
De- Therefore, the electron-electron scattering rate in- 
creases much faster than the electron-phonon scattering 
rate, which makes the electron system thermal and the 
TTM description valid. ^5- According to Kaganov et a/.'s 
calculation, 12 ^ gQgg ^s g{T) ~ at low 

T and g{T) constant at high T. For normal metals, 
since Ce ^ Cp even down to 10 K, the TTM relation 
Eq. ^ gives r ~ T~^. In the case of heavy fermion 
compounds, Ce is comparable to or larger than Cp at 
low T, and therefore Cp ^ for the Debye phonons 
produces ~ T as the leading T dependence. Indeed, 
the low-T divergence of the calculated t in Fig. |^ can 
be well described as T~^. Therefore, the discrepancy 
with experimental data cannot be fixed by simply chang- 
ing parameters, and it indicates that there is additional 
physics involved other than the large DOS at the Fermi 
level to explain a 100-fold increase of r below Tk- 

If we take the approximation that the difference be- 
tween Ep and the bottom of the heavy electron band is 
about Tk, then from the isotropic parabolic heavy elec- 
tron dispersion model and the carrier density of n = 0.85 
measured by the Hall effect for YbAgCu4,l'' we can es- 
timate Vp A km/sec. Ultrasonic measurements show 
that the longitudinal sound velocity along the [111] direc- 
tion is about 4.4 km/sec This comparison shows that 
in fact the two velocities are comparable to each other, 
supporting the possibility of vp < Vs- We should note 
that experiments like de-Haas van-Alphen are required 
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for a more reliable estimate oi vp. Therefore, we hy- 
pothesize that the Fermi velocity is lower than the sound 
velocity, and discuss how this affects the electron-phonon 
relaxation dynamics. 

As depicted in Fig. ^1 if the Fermi velocity is slower 
than the sound velocity [Fig. Iior b)]. then momentum 
and energy conservation requirements prohibit scatter- 
ing between heavy electrons and phonons. In the actual 
heavy fermion DOS, heavy electron dispersion exists only 
within the peak DOS around Ep, and electrons outside 
the peak have a regular light mass. Therefore, the above 
blocking mechanism applies only when both the initial 
and final electron states are within the peak of the DOS. 
If one or both of the states are outside the peak, regu- 
lar electron-phonon scattering is expected. These effects 
can be simply incorporated into the simulation by block- 
ing electron-phonon scattering within the peak DOS. The 
results of the simulation are shown in Fig. llll as solid cir- 
cles, for which the blocked energy range is from -24 meV 
to 24 meV, K^p = 0.23 eV/ps, and Kee = 700 x K^p. 
This Kep value is smaller than the value for LuAgCu4 
case, by about half in terms of the scattering matrix ele- 
ment. It shows reasonable agreement with experimental 
data (open circles). 

Since the electron system has a thermal distribution, 
the TTM approach, with a similar blocking of electron- 
phonon scattering within the peak, can be used to de- 
scribe the relaxation between the electrons and phonons, 
which was discussed in Ref. 0. In this TTM approach, 
we also included the possibility that some electron- 
phonon scattering within the peak is allowed due to the 
anisotropy of vp{k) and Vs{k), that is, vp > Vs along 
some directions and vp < Vg along other directions. The 
results were shown as a solid line in Fig. 3 in Ref. 
and reproduced in Fig. ^] (solid line), which shows good 
agreement with experimental datai^ 

The physical idea behind the increase of r below Tr- 
is simple: As T decreases below Tr-, or equivalently 
the width of the DOS peak, since the electron system 
has a thermal distribution, the main electron-phonon 
relaxation processes should occur within the peak [see 
Fig. inia)] , where the momentum and energy conserva- 
tion laws block the electron-phonon scattering. There- 
fore, the electron-phonon relaxation time increases very 
rapidly as T is lowered below The results show 

that complete or substantial blocking of electron-phonon 
scattering processes within the DOS peak is essential to 
explain the rapid increase of the electron-phonon relax- 
ation time below Tk in YbAgCu4. 

IV. SUMMARY 

We provide a theoretical analysis of the ultrafast 
relaxation dynamics observed by femtosecond time- 
resolved optical spectroscopy in isostructural LuAgCu4 
and YbAgCu4; the former is a normal metal and the lat- 



ter is a hea vy fermion compound. For norma l metals, we 
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FIG. 10: Schematic pictures explaining different electron- 
phonon scattering processes for (a) vp > Vs and (b) vf < Vs 
cases. If wf < "Us, then the electron dispersion changes more 
slowly than the phonon dispersion. Therefore, the distance 
between the two fc-space spheres with energies e and e — u is 
larger than the phonon momentum q, and no momentum and 
energy conserving scattering process is possible, as shown in 
(b) in this figure. 



find that a nonthermal electron distribution is responsi- 
ble for a temperature-independent electron-phonon relax- 
ation time at low temperatures, instead of a diver- 
gent behavior predicted by the two-temperature model. 
For heavy fermion compounds, we find that prohibiting 
electron-phonon scattering within the density-of-states 
peak near the Fermi energy is crucial to explain the rapid 
increase of the electron-phonon relaxation time below the 
Kondo temperature. On the basis of the estimated Fermi 
velocity and the measured sound velocity, we propose the 
hypothesis that the slower Fermi velocity compared to 
the sound velocity provides this blocking mechanism due 
to energy and momentum conservation laws. We find 
good agreement between the experimental data and our 
model for both normal metals and heavy fermion com- 
pounds. 

We thank V. V. Kabanov for useful discussions. This 
work has been supported by U.S. DOE. 
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FIG. 11: Solid circles: relaxation time calculated with peaked 
DOS and blocking of e-p scattering within the peak DOS due 
to vf < Va. Open circles: experimental data for YbAgCu4. 
Solid line: results obtained in Ref. using TTM with partial 
blocking of the e-p scattering within the peak. 
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